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Abstract 



We perform a smoothed analysis of the GCC-condition number 'if (A) 
of the linear programming feasibility problem 3x £ E m+1 Ax < 0. 
Suppose that A is any matrix with rows <Zj of euclidean norm 1 and, 
independently for all i, let Oj be a random perturbation of 5j fol- 

■ lowing the uniform distribution in the spherical disk in S m of angu- 
k*" | lar radius arcsiner and centered at a,. We prove that EQn^A)) = 

0(mn/a). A similar result was shown for Renegar's condition num- 
ber and Gaussian perturbations by Dunagan, Spielman, and Teng 

■ [arXiv:cs.DS/0302011]. Our result is robust in the sense that it easily 

extends to radially symmetric probability distributions supported on 

■ a spherical disk of radius arcsiner, whose density may even have a sin- 
| gularity at the center of the perturbation. Our proofs combine ideas 

■ from a recent paper of Biirgisser, Cucker, and Lotz (Math. Comp. 77, 

No. 263, 2008) with techniques of Dunagan et al. 

i . 
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1 Introduction 

A distinctive feature of the computations considered in numerical analysis is 
that they are affected by errors. A main character in the understanding of 
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the effects of these errors is the condition number of the input at hand. This 
is a positive number measuring the sensitivity of the output with respect 
to small perturbations of the input. The best known condition number is 
that for matrix inversion and linear equation solving, which takes the form 
k(A) = \\A\\ WA^ 1 1| for a square matrix A. Condition numbers not only occur 
in round-off analysis, but also appear as a parameter in complexity bounds 
for a variety of iterative algorithms in linear algebra, linear and convex 
optimization, and polynomial equation solving. Yet, condition numbers are 
not easily computable. As a way out for this situation, Smale suggested to 
assume a probability measure on the set of data and to study the condition 
number of this data as a random variable. Examples of such results abound 
for a variety of condition numbers. For more details and references we refer 
to Smale's survey [25] and the recent survey [4]. 

Renegar [20, 21, 22] was the first to realize that the computational com- 
plexity of linear programming problems can be bounded by a polynomial in 
the number of variables and inequalities and a certain condition measure of 
the input. This condition measure is given by the inverse distance of the 
input to the set of ill-posed systems. More specifically, it is well known that 
for a given matrix A £ R riX ( m+1 ), either the system Ax < or its dual sys- 
tem A T y = 0, y > have a solution, unless we are in an ill-posed situation. 
The (homogeneous) linear programming feasibility problem is to decide this 
alternative for given A and to compute a solution of the corresponding sys- 
tem. A primal-dual interior point method is used in [12] to solve the linear 
programming feasibility problem within 

0(Vm + n (ln(m + n) + ln^(A))) (1) 

iterations, with each step costing at most 0((m+n) 3 ) arithmetic operations. 
Here, the GCC- condition number ^{A) is a variant of Renegar's condition 
number introduced by Goffin [17], and later generalized by Cucker and Che- 
ung [8] (see §2.3 for the definition). The advantage of ^{A) over Renegar's 
condition number is that this quantity can be neatly characterized in terms 
of spherical geometry, which greatly facilitates its probabilistic analysis. 

Thus the running time of the primal-dual interior point method used 
in [12] is controlled by 'tf(A). Understanding the average-case behaviour 
of this algorithm therefore boils down to studying In'rf(A) for random in- 
puts A. Motivated by this observation, a lot of efforts have been devoted to 
the average-case analysis of the random variable ^(A), i.e., to compute the 
expected value (or the distribution function) of In ^(A) for random matri- 
ces A. In most cases, the matrices A are assumed to have random entries 
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which are i.i.d. standard normal. We remark that since the condition number 
^(A) is multi-homogeneous in the rows dj of A, this is equivalent to con- 
sidering ^(A) in the case where a±, . . . , a n are i.i.d. uniformly distributed in 
unit sphere S m := {x G R m+1 \ \\x\\ = 1}. 

The papers dealing with the average analysis of ft (A) are easily sum- 
marized. A bound for E(ln^(>l)) of the form 0(min{n,mlnn}) was shown 
in [9]. This bound was improved in [13] to maxjln m, In In n} + 0(1) as- 
suming that n is moderately larger than m. Still, in [10], the asymptotic 
behavior of both ft (A) and ln^(A) was exhaustively studied and these re- 
sults were extended in [18] to matrices A G {S m ) n drawn from distributions 
more general than the uniform. Finally, in [5], the exact distribution of 
^(A) conditioned to A being feasible was found and asymptotically sharp 
tail bounds for the infeasible case were given. In particular, it was shown 
that E(ln c ff(A)) = O(lnm). Our method yields another proof of this re- 
sult (Cor. 1.3). 

1.1 Smoothed analysis 

The problem of average-case analysis is that its results strongly depend on 
the distribution of the inputs, which is unknown, and usually assumed to 
be Gaussian for rendering the mathematical analysis feasible. Spielman 
and Teng [26, 27, 28] suggested in 2001 the concept of smoothed analysis 
of algorithms, which is a new form of analysis of algorithms that arguably 
blends the best of both worst-case and average-case. They used this new 
framework to give a more compelling explanation of the speed of the simplex 
method (for the shadow- vertex pivot rule). 

The general idea of smoothed analysis is easy to explain. Let T : W — > 
IR + U{oo} be any function (measuring running time, log of condition numbers 
etc.). Instead of showing "it is unlikely that T(a) will be large," one shows 
that "for all a G W and all slight random perturbations a of a, it is unlikely 
that T(a) will be large." If we assume that a is multivariate normal with 
mean a and variance a 2 , in short a G N(a,a 2 ), then the goal of a smoothed 
analysis of T is to give good estimates of 

sup Prob {T(a) > e -1 }. 

aGRP a€N(a,a 2 ) 

In a first approach one may focus on bounds on the expectations. 

For many situations of interest, it turns out that in smoothed analysis, 
there is only a weak dependence on the chosen model of random perturba- 
tions. A first formal evidence of this robustness phenomenon was given by 



3 



Cucker, Hauser, and Lotz [11] for certain conic condition numbers, stated 
below as Theorem 1.1. 

Dunagan et al. [15] (see also Spielman and Teng [27]) performed a 
smoothed analysis of the running time of interior point methods of linear 
programming by analyzing Renegar's condition number c rf / , a variant of the 
GCC condition number % '. Among other things, they proved the following: 

sup E AeN{A ^ } (ln^(A)) = O(ln^). 

Here the supremum is over all A G M nx ( m+1 ) of Frobenius norm at most 
one. Our main result (Theorem 1.2) shows that a similar bound actually 
holds for a large class of random perturbation laws. 



1.2 A geometric approach to conic condition numbers 

Biirgisser et al. [6, 7] recently extended a result of Demmel [14] on conic 
condition numbers from average-case analysis to smoothed analysis. There, 
the perturbations of the inputs are modelled by uniform instead of Gaussian 
distributions. This allows to perform the analysis in a general geometric 
framework that we explain next. 

The set of ill-posed inputs to a computational problem is modelled as 
a lower dimensional subset E of the data space @, which is assumed to be 
furnished with a metric d, a distance function dist (not necessarily a metric), 
and a volume measure. In our cases of interest, £F is a Riemannian man- 
ifold, d is the corresponding metric, and dist = sin(d). The corresponding 
condition number ^(a) of an input a G £F is then defined as 

tf(o) = - } — . (2) 
v ; dist(a,E) v ' 

This is an appropriate definition for many applications. In this model, the 
set of inputs a with condition ^(a) > e _1 is given by the e- neighborhood 

T dist (E, £ ) = {a G 9 | dist(a,S) < e}. 

Let B(a,a) := {a G Q> \ d(a,a) < a} denote the ball centered at a G 3> of 
radius a. The task of a uniform smoothed analysis of consists of providing 
good upper bounds on 

sup Prob {'if (a) > e -1 }, 

a£@a<=B(a,a) 
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where a is assumed to be chosen uniformly at random in B(a,a). The 
probability occurring here thus has an immediate geometric meaning: 




ob ma)>e- 1 } = 

a,a) 



vol (r dist (£,£)n£(a,a)) 
vol (B(a,a)) 



(3) 



Thus uniform smoothed analysis means to provide bounds on the volume 
of the intersection of e- neighborhoods of S, relative to the distance func- 
tion dist, with balls of radius a. We note that for compact data spaces 3$, 
uniform smoothed analysis interpolates transparently between worst-case 
and average-case analysis. Indeed, when a = we get worst-case analysis, 
while for a = diam(^) we obtain average-case analysis. 

For conic condition numbers, this general concept specializes in the fol- 
lowing way. The data space is the unit m-sphere S m C R m+1 , £ is a lower 
dimensional subset of S m such that £ = — £ (in many applications it is an 
algebraic hypersurface) and the conic condition number ^(a) of a £ S m is 
defined by (2) with respect to the distance function dist(a,fe) := sin d(a, b), 
where d refers to the angular (i.e., Riemannian) distance in S m . (We could 
as well consider the data space as the real projective space P m on which dist 
defines a metric.) 

For defining the GCC-condition number, one takes as the data space the 
nth power (S m ) n = S m x • • • x S m of the sphere S m with the metric defined as 
d(A,B) := maxi<j< n ii(ai,6i), where A = (a 1 , . . . ,a n ), B = (bi, . . . , b n ) and 
d(cii,bi) again denotes the angular distance in S m . The distance function is 
defined by d\st(A,B) := sind(A, B) and the GCC-condition number ^(A) 
is given by ^(A) = l/dist(A, S n , m ), where Y> n , m C (S m ) n denotes the set of 
ill-posed inputs, which is a semialgebraic subset of codimension one (cf. §2.3 
for details). 

1.3 Adversarial distributions 

An advantage of the uniform model is that, for conic condition numbers, 
results on smoothed analysis easily extend to more general families of prob- 
ability distributions. To obtain such robustness results we rely on a general 
boosting technique developed in Hauser and Miiller [18] and apply it simi- 
larly as in Cucker et al. [11] 

The framework is the following (cf. §2.4). Consider the spherical cap 
B(a, a) := {x £ S m \ d(x, a) < a} in the sphere S m centered at a £ S m and 
having angular radius a G (0,7r/2]. Let v denote the uniform measure on 
B(a, a) and suppose that [ig- is a z/- absolutely continuous probability measure 
on S m , say /ia(G) = f G fdv for measurable sets G. We further assume that 
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the density / is of the form f(x) = g(smd(x,a)) with a monotonically 
decreasing function g: [0,a] — > [0,oo] of the form 

g (r) = Cr~ fi h(r), 

where a := sin a, (3 < m, h : [0, a] — > [0, oo) is a continuous function satis- 
fying h(0) 7^ 0, and C = I m (a)/I m -p(a) is a normalizing constant, where 
Ifc(a) := J a (sint) fc_1 dt, cf. §2.2. Thus, the support of /x„ is contained in 
B(a,a) and the density of fi^ is radially symmetric with a pole of order f3 
at a. Such distributions were called adversarial in [11]. The exponent f3 
and the quantity H := sup 0<r<CT h(r) are the only parameters entering the 
bound below. In the special case = 0, the density of does not have a 
singularity and the situation considerably simplifies: We have C = 1, g = h 
and the estimate Ha(B) < H ■ v{B) holds for any measurable set G. 

To extend this estimate to the general case, the smoothness parameter s 
of a ^-absolutely continuous distributions was defined in [18], and in [11, 
Lemma 3.2], it was shown that s = 1 — (3/m. This means that s is the 
largest number s' > for which the following is true: For every e > 
there exists 5(e) > such that v(G) < 5(e) implies Ha(G) < u(B) s ~ e for 
all measurable sets G. This allows to obtain tail bounds for /x„ from tail 
bounds for v. 

Biirgisser et al. [7] provided a general smoothed analysis of conic con- 
dition numbers for uniform perturbations. This was recently extended by 
Cucker et al. [11] to the model of adversarial perturbations, who obtained 
the following robust smoothed analysis estimate. 

Theorem 1.1 Let & be a conic condition number with set of ill-posed 
inputs S C S m and assume that E is contained in an algebraic hypersurface 
of degree d. Then we have for a random perturbation from any adversarial 
distribution fj,a on B(a,a) with center a £ S m , angular radius a £ (0, ir/2] 
and parameters f3,H,a = sin a that 

m /i r~s>\ , rn 2 d 1 , 2eH 2 m , 13vr 

E„_ n^f < In + -—In— — + In . 

^ aV ' ~ a 1-p/m ln(7rm/2) 2 

It is remarkable that the only problem dependent parameters entering 
the above bound are the dimension m and the degree d. The only distribu- 
tion dependent parameters entering the bound are a, /3, and H. This result 
has a wide range of applications to linear and polynomial equation solving. 
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1.4 Main results 



The goal of this paper is to prove the following analogue of Theorem 1.1 for 
the GCC-condition number of the linear programming feasibility problem. 

Theorem 1.2 Let ^ denote the GCC condition number defined on (S m ) n 
and n > m + 1. Suppose that di is randomly chosen from an adversarial 
distribution u.^ on Bfa, a) with center ai £ S rn , angular radius a 6 (0, 7r/2], 
and parameters f3,H,cr = sin a, independently for i = 1, . . . , n. Then the 
random matrix A with rows satisfies 



In the case f3 = 0, we have more precisely, 

E(ln^(A)) < 121nn+171nm + 61ni + 81n# + 29 = 0(ln^). 

As an application let T{A) denote the number of iterations of the primal- 
dual interior point method of [12] for solving the linear programming fea- 
sibility problem 3x 6 M m+1 Ax < (or its dual problem). Theorem 1.2 
implies that 



for a random matrix A with independent rows a« from a adversarial distri- 
bution fiai with parameters f3, H. 

For the uniform distribution on S m , by essentially the same method, 
we can improve the estimates of Theorem 1.2 obtaining a result that was 
previously shown in [7] by a very different technique. 

Corollary 1.3 Suppose that the rows of the matrix A are independently 
chosen in S m according to the uniform distribution and n > m + 1. Then 



The paper is organized as follows. Section 2 is devoted to preliminaries on 
spherically convex sets, the GCC-condition number, and adversarial proba- 
bility distributions. In Section 3 we state and prove probability tail bounds 
for the GCC-condition number under adversarial random perturbations and 
prove our main results. The proof essentially reduces the adversarial to the 
uniform case and then uses geometric arguments for the uniform case. A 
principal ingredient of the proof of the uniform case is an upper bound on 




E(T(A)) = 0((l-iy 2 .Vn-.ln^L), 



E(ln^(A)) = O(lnm). 
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the volume of the neighborhood of spherically convex sets (Theorem 3.3), 
that is stated in Section 3.1, but whose proof is deferred to Section 4. The 
proof of the latter result proceeds along the lines of [7] and uses some deeper 
results from integral and differential geometry (Weyl's tube formula, kine- 
matic formula). 

Acknowledgments. We thank Felipe Cucker and Martin Lotz for nu- 
merous helpful discussions. We are grateful to an anonymous referee for 
constructive criticism that led to more general results and a significantly 
better presentation of the paper. 

2 Preliminaries 

2.1 Convex sets in spheres 

A general reference about convex sets is [30]. Glasauer's thesis [16] is a 
useful reference for the integral geometry of spherically convex sets. 

A convex cone in M m+1 is a subset that is closed under addition and 
multiplication with nonnegative scalars. We denote by cone(M) the convex 
cone generated by a subset M C ]R m+1 . More specifically, the convex cone 
generated by points a±, . . . , a k G E m+1 is given by 

cone{ai, . . . , a k } := {x £ R m+1 \ 3X 1 > 0, . . . , X k > x = Y%=i A * a 0- 

A convex cone C is called pointed iff C D (— C) = {0}. It is known that 
C is pointed iff C \ {0} is contained in an open halfspace whose bounding 
hyperplane goes through the origin. Clearly, if a« 7^ for i = l,...,k, 
then conejai, . . . ,a&} is pointed iff is not contained in the convex hull 
conv{ai, . . . ,a k }. 

We use convex cones to define the notion of convexity for subsets of the 
sphere S m = {x G W n+1 \ \\x\\ = 1}. Let x,y G S m be such that x / ±y. 
We call [x, y] := conejx, y}nS m the great circle segment connecting x and y. 

Definition 2.1 A subset K of S m is called (spherically) convex iff we have 
[x, y] C K for all x,y £ K with x 7^ ±y. We call K properly convex iff it is 
nonempty, convex, and does not contain a pair of antipodal points. 

We denote by sconv(M) := cone(M) n S m the convex hull of a subset M 
of S m , which is the smallest convex set containing M. Clearly, M is convex 
iff M = sconv(M). The closure of a convex set is convex as well. It is easy to 
see that a convex subset K of S m is contained in a closed halfspace, unless 
K = S m . A properly convex set K is always contained in an open halfspace. 
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Definition 2.2 The dual set of a subset M C S m is defined as 
M := {a G S m | Vx G M (a, x) < 0}. 

Clearly, M is a closed convex set disjoint to M. The hyperplane sepa- 
ration theorem implies that the dual of M equals the closure of sconv(M). 
We note that M C N implies M D N. Finally, it is important that M has 
nonempty interior iff M does not contain a pair of antipodal points, that is, 
"nonempty interior" and "properly convex" are dual properties. 

By a convex body K in S m we will understand a closed convex set K such 
that both K and K have nonempty interior, i.e., both are properly convex. 
The map K i-> K is an involution of the set of convex bodies in S m . 

2.2 Distances, neighborhoods, and volumes 

We denote by d(a, 6) G [0, ir] the angular distance between two points a, b 
on the sphere S m . Clearly, this defines a metric on S m . The (closed) ball of 
radius a G [0, n] around a G S m is defined as 

B(a, a) := {x G S m \ d(x, a) < a} = {x G S m \ (a, x) > cos a}. 

This is the same as the spherical cap with center a and angular radius a. 
B(a, a) is convex iff a < ir/2 or a = n. In the case a < n/2, the dual set of 
B(a, a) equals B{— a, it/2 — a). 

For a nonempty subset M of S m we define the distance of a G S m to M 
as d(a, M) := inf{d(a, x) \ x G M}. The dual set of M can be characterized 
in terms of distances by: a G M <^=^ d(a, M) > tt/2. 

Lemma 2.3 Let K he a convex body in and a G S m \ (K U K). Then 
d(a, K) + d(a, K) = tt/2. 

Proof. Let b £ K such that 99 := d(a,b) = d(a,K). Since a G" if we have 
99 < 7r/2. The point 6* := (a, 6)6 is therefore nonzero and contained in 
C := cone(K). Put p* ■= a - b*. Then (p*,6) = 0, (p*,o) = sin 2 ip, and 
{p*,p*) = sin 2 99. In particular p* 7^ 0. 

By construction, 6* is the point of C closest to a. It follows that {x G 
fljm+i J (jf^x^ = 0} is a supporting hyperplane of C. Hence (p*,x) < for 
all x G C and the point p := p*/\\p*\\ therefore belongs to K. Moreover, 
(p, a) = simp, which implies d(a,p) = tt/2 — ip. Hence 

d(a, K) + d(a, K) < d(a, 6) + d(a,p) = tt/2. 
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To complete the proof it suffices to show that d(a,K) = d(a,p). Suppose 
there exists p' G K such that d(a,p') < d(a,p). Then d(b,p') < d(b,a) + 
d(a,p') < d(b, a) + d(a,p) = tt/2 which contradicts the fact that b G K. □ 

Sometimes it will be useful to consider the projective distance between 
two points a, b G S m , which is defined as dp(a, b) := sind(a, b). It is straight- 
forward to check that dp satisfies the triangle inequality. However, it is not 
a metric on S m , as d$>(a,b) = iff a = ±b. Hence the ball of radius sin a, 
measured with respect to the projective distance, equals B(a, a)L) B (— a, a) . 
We denote this set suggestively by B(±a,a) and call it the projective ball 
with center ±a and radius a. 

For < ip < tt/2, the ip -neighborhood of a nonempty subset M of S m is 
defined as T(M, ip) := {x G S m \ d(x,M) < ip}. If M is the boundary dK 
of a properly convex set K in S m , we call 

T (dK, ip) := T(dK, <p)\K and T^dK, <p) := T(dK, ip) n K 

the outer ip -neighborhood and inner ip -neighborhood of clfT, respectively. 
Clearly, we have T(dK, ip) = T Q (dK, ip) U T-{dK, ip). 

In order to compute the m-dimensional volume of such neighborhoods, 
the following functions J m ,k(ct) are relevant: 

J m ,k( a ) ■= /"(sinp)^ 1 (cos p) m - k dp (1 < k < m). (4) 
J o 

Recall that O m := vo\S m = 2vr( m+1 )/ 2 /r((m+l) /2) equals the m-dimensional 
volume of S™. It is known that vo\T(S m ~ k ,a) = O m -kOk-iJm,k{o{). Some 
estimations of these volumes can be found in [7, Lemmas 2.1-2.2]. 

2.3 The GCC condition number 

We study the problem of deciding for a given instance A G IR rax<m+1 ) whether 
there exists a nonzero solution x G R m+1 \ {0} such that Ax < 0. In the 
following we assume that n > m + 1. (The other case is considerably less 
interesting.) Without loss of generality we may assume that the row vectors 
a,i have euclidean length one, and hence interpret A = (ai,...,a n ) as an 
element of the product (S m ) n of spheres. 

We write sconv(A) := sconvjai, . . . , a n } for the convex hull of the given 
points. The set of solutions in S m of the system of inequalities Ax < 
equals the dual set of sconv( J 4). 
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Definition 2.4 An instance A G (S' m ) n is called feasible iff its set of solu- 
tions is nonempty, otherwise A is called infeasible. An instance A is called 
strictly feasible iff its set of solutions has nonempty interior. We denote by 
T n<m and m the set of feasible and strictly feasible instances, respectively. 
The set of ill-posed instances is defined as S„ jm := T n ^ m \ T^ m . The set of 
infeasible instances is denoted by I n ,m- 

Remark 2.5 An instance A G (S m ) n is strictly feasible iff sconv(A) is prop- 
erly convex, that is, cone(A) is pointed. Furthermore, a feasible instance A 
is ill-posed iff is contained in the euclidean convex hull of a±, . . . , a n (cf. [5, 
Lemma 3.2]). 

We remark that J- n ,m is a compact subset of (S m ) n with nonempty inte- 
rior T° m and topological boundary £ nim . Moreover, I n ,m is nonempty and 
S„ jm is also the topological boundary of X n , m . (Here we use n > m + 1.) 

We define a metric on (g m ) n by setting for A,Be (S m ) n with compo- 
nents aj, bi G S m 

d(A, B) := max d(ai,bi) . 

l<i<n 

The distance of A to a nonempty subset M C (5 m ) n is defined as cZ(^4, M) := 
mf{d(A, B)\B G M}. We denote by 5(1, a) := {^4 G (S™) n | d(A, A) < a} 
the closed ball with center A and radius a. Clearly, this is the product of 
the balls B(a>i, a) for i = 1, . . . , n. 

The following definition is due to Goffin [17] and Cheung and Cucker [8]. 

Definition 2.6 The GCC condition number of A G (S m ) n is defined as 
tf(A) = l/sind(A£ n , ro ). 

This condition number can be characterized in a more explicit way. 

Definition 2.7 A smallest including cap (SIC) for A G (£"")" is a spherical 
cap of minimal radius containing the points a±, . . . , a n . 

We remark that by a compactness argument, a SIC always exists. It can 
be shown that a SIC is unique if A is strictly feasible. However, for infea- 
sible A, there may be several SICs (consider for instance three equidistant 
points on the circle). We denote the radius of a SIC of A by p(A). An 
instance A is strictly feasible iff p(A) < tt/2. For more information on this 
we refer to [10, 5]. 

The following result is due to Cheung and Cucker [8] . This characteriza- 
tion is essential for any probabilistic analysis of the GCC condition number. 
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Theorem 2.8 We have 



d(A, X n ,m) — 



§ - p(A) if A G T n , m , 
p(A)-\ ifAe(S m r\T n , 



In particular, d(A, < | and ^(A) 1 = | cos p(A)\. 

2.4 Adversarial probability distributions 

The proof of our robustness result relies on a general boosting technique 
developed in [18], that allows to extend probability tail bounds obtained 
for one fixed distribution to larger classes of distributions. We explain this 
technique in our situation of interest. 

Let v := z/a jCr denote the uniform distribution on the spherical disk B(a, a), 
where a G (0, §] and a := sin a. We assume that p is a zv-absolutely con- 
tinuous probability measure, i.e., it can be written with a density / as 
p(G) = J G fdv for Borel measurable sets G. In certain cases, it is possible 
to bound p(G) in terms of v(G) if the latter is sufficiently small. This is done 
with the smoothness parameter s of p with respect to v, which is defined 
as s := lim^o inf 5, where we have set for 5 G (0, 1), using the convention 
InO := — oo, 



If s is positive, we say that p is uniformly v- absolutely continuous. In this 
case, it is easy to see that s is the largest nonnegative real number s' with 
the property that for all e > there exists a tolerance 5(e) > such that 
v(G) < 5(e) implies p(G) < v(G) s '~ £ for all G, cf. [18]. 

We will apply this framework to a specific class of distributions p. Not 
only will it be important to know the smoothness parameter, but also to 
explicitly compute bounds for the tolerance 5(e). 

An adversarial probability distribution p^, for a G S m , was defined in [11] 
as a ^-absolutely continuous measure given by Pa(G) = J G fd u, where again 
v denotes the uniform distribution on B(a,a). We further require that the 
density / is of the form f(x) = g (sin d(x, a)) with a monotonically decreasing 
function g: [0, a] — > [0, oo] given by 



where < /3 < m, and h: [0, a] — > [0, oo) is a continuous function satisfying 
h(0) 7^ 0. Thus, the support of p-a is contained in B(a,a) and the density 
of pa is radially symmetric with a pole of order (3 at a. Clearly, a is the 




g(r) = Or-? h(r), 
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mean of a with respect to jia- It is convenient to assume the normalization 
C := I m (a)/I m ^fs(a), where 

h{a) := J k ,k( a ) = / (sin i)^ 1 (it, 
J o 

compare (4) (we slightly deviate here from the notation in [11]). Then /i„ is 
a probability distribution iff 

pa 

I m -p(a)= (sin t) m ~ ^ h{smt) dt. 
Jo 

The maximum of h is denoted by H := sup 0<r<cr h(r), which is easily seen to 
satisfy H > 1. We note that the uniform distribution on B(a, a) is obtained 
by choosing /3 = 0, C = 1 and for g = h the function identically equal to 1. 

The following technical lemma is an immediate consequence of [1 1 , Lem- 
mas 3.2-3.3, Equation (3.1)]. 

Lemma 2.9 1. The smoothness parameter of \i with respect to v equals 
s = 1 — f3/m. 

2. For c := = \(\ — ^) < \ and the tolerance 5 C defined by 




we have for all measurable G C S m , 

u(G) < 8 C =► Ha{G) < (u(G)f. 

Remark 2.10 In the case = where the density of fia has no singularity 
the situation simplifies. Clearly, /ia(G) < Hv{G) for all G. This directly 
implies that the smoothness parameter equals 1. Moreover, the tolerance 
5 = is sufficient for the implication v{G) < 5 =>■ ^(G) < z/(G)2. 

Since we will work in the product of spheres (S m ) n we define on it 
the adversarial distributions fi^ with center A = (ai,...,a n ) G (S m ) n by 
taking the product measure fi^ := /Ltai x ... x ^a n . This is a probability 
distribution whose support is contained in the product of caps B(A,a) : = 
B(a±, a) x ... x B(a n , a). 
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3 Robust smoothed analysis of ^(A) 



The goal of this section is to provide smoothed analysis estimates for the 
condition number ^(A) in the model where A = (a±, . . . ,a n ) is chosen at 
random according to an adversarial distribution fi^. The center A G (S m ) n 
of the perturbation is arbitrary. Recall that F n ,m an d ^n,m denote the sets 
of feasible and infeasible instances A 6 (S m ) n , respectively. 

Theorem 3.1 Let a G (0,7r/2], a = sin a, and A G (S m ) n , wheren > m + 1. 
Assume that A G -B(j4,a) is chosen at random according to an adversarial 
distribution /x^. Let the exponent c = ^(1 — fi/m) and the tolerance 5 C be 
defined as in Lemma 2.9. Then we have 

Prob{A G T n , m and V(A) > t} < n ^ 13m (™ + V ^j t ~c (F) 

provided t > 13 ™^ +1 - ) . Moreover, we have for t > 1, 

^ , r . ^ , . , ( 1690m 2 (m + 1) \ c . . . . 
Prob{A G l n , m and <g(A) > t} < n t~ c (S~ c + mint). 



V 4a 2 

(I) 

Remark 3.2 For the uniform distribution (/3 = 0, h = 1) the bounds in (F) 
and (I) can easily be improved by avoiding the use of Lemma 2.9. In partic- 
ular, on the right-hand side of (F) and (I) one gets t^ 1 and t^ 1 hit instead 
oft" 1 / 2 and r^lnt 

The overall strategy of the proof of Theorem 3.1 is the same as in [15]. 
However, the crucial component in [15], namely a result due to Ball [2], is 
substituted by Corollary 3.4, which is stated in the next section. 



3.1 A bound on the volume of convex neighborhoods 

We state here an upper bound on the volume of the intersection of neigh- 
borhoods of spherically convex sets with spherical caps. 

Theorem 3.3 Let K be a properly convex subset of S m , let a G S rn , and 
< a,(p < 7r/2. Then, writing a = sin a and e = sin^, we have the 
following upper bound for the volume of the outer neighborhood of dK: 

vo\(T (dK, <p)r\B(±a,q)) < y fm\ A £\ m ~* ( £ \ k , m °™ /£\ m 
vo\B(±a,a) ~ ^ W ^ °' ^ 20 m _i W 
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The same upper bound holds for the volume vol( - T ' ^^± c ^ za ' a ^ of the inner 
neighborhood of dK. 

The main result of [7, Theorem 1.2] gives a bound on the volume of the 
intersection of a projective ball B(±a, a) = B(a, a) U B(—a, a) with the <p- 
neighborhood of a real algebraic hypersurface in the sphere S m , given as the 
zero set of polynomial of degree d. The above theorem says that essentially 
the same volume bound holds for the boundary of a properly convex set K 
in S m , if we formally replace in this bound the degree d by 1/2. 

The proof of Theorem 3.3, which is quite involved, is deferred to Sec- 
tion 4. By essentially the same argument as in the proof of [7, Prop. 3.5] 
one can derive from Theorem 3.3 the following corollary. 

Corollary 3.4 Under the assumptions of Theorem 3.3 we have the follow- 
ing upper bound for the volume of the outer neighborhood of dK: 

vo\(T (dK,ip)nB(±a,q)) < 13m e_ < ^_ 

vo\B(±a,a) ~ 4 a £ ~ 2m ' 

The same upper bound holds for the relative volume of the inner neighbor- 
hood of dK. 



3.2 Two auxiliary results 

The proofs of the following two results are similar as in Dunagan et al. [15]. 
We use the notation [n] := {1,2,..., n} for n G N. 

Proposition 3.5 Let A= (ai, . . . , a n ) G J 7 " m , < ip < tt/2, and e = sin ip. 
If^(A) > (m + 1) e~ x , then there exists i G [n\ such that 

ai e T{dKi,p)\Ki, 

where K { := -sconv{ai, . . . , aj_i, a i+ i, ... ,a n }. 

Proof. There exists q G sconv(A) such that (ai,q) > for all i G [n\. 
Indeed, if q is taken as the center of the SIC of A then this follows from [10, 
Lemma 4.5] (see also [5, Lemma 3.2]). 

We note that Oj K; L for all i G [n]. Otherwise G convjai, . . . ,a n }, 
hence A G X ra ,m> which contradicts our assumption that A is strictly feasible. 
It follows that d(di,dKi) = d(cii,Ki). 

We assume now d(ai,Ki) > p for all i G [n]. Our goal is to show that 
sind(A, S n>m ) > ^pye- Then we are done, since ^{A)^ 1 = sin d(A, S n;m ) 
by Definition 2.6. 
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We proceed now similarly as in [15, Lemma 2.3.10]. By continuity we 
assume w.l.o.g. that ip < tt/2. We distinguish two cases. If G" Ki, then 
Lemma 2.3 tells us that d{a,i,Ki) +d(ai,Ki) = tt/2. Hence d(ai,Ki) < 
tt/2 — if. Choose pi G int(ifj) such that d(a,i,pi) < tt/2 — if. This implies 
(cii,pi) > cos(-7r/2 — ip) = e. In the case etj G Ki we take any pi G int(iiQ) 
close enough to such that (ai,pi) > e. 

In both cases we have achieved the following 

(ai,pi) > e and Vj ^ i (aj,pi)>0. (5) 

This implies for all i that (pi, q) > 0, as q G conejai, . . . , a n }. 
Consider now for i G [n] the following convex sets in S m 

Q := {x e S m \ (di, x) > — - — and (x, q) > 0} 

m + 1 

containing p^. We claim that the intersection of any m + 1 of these sets 
is nonempty. Indeed, let / C [n] be of cardinality m + 1 and consider 
P* '■= J2jeiPj- Note that ||p*|| < 1. We obtain for any i G /, using (5), 

= — xry^(oi>Pj) > — l —r{ai,Vi) > £ , ; • 
m + l^ m + 1 m + 1 

Moreover, (p*, q) > 0, hence p* / 0. It follows that p := is contained 

in Ci for any i £ I, which shows the claim. 

Consider the affine hyperplane E := {x G R m+1 | (x,q) = 1} of dimen- 
sion m and the perspective map 

tt : {x G S m | (x, g) > 0} -+ E 1 , x i-)- (5, x) _1 x. 

Then the 7r(Ci) are convex subsets of E, with the property that any m+1 
of these have a nonempty intersection. Helly's theorem [30] implies that 
tt{C\) fl • • • fl ir(C n ) is nonempty. Hence there is a point a G D^Li We 
have d(a,i,a) < a := arccos((m + l) _1 e) for all i G [n]. Hence the spherical 
cap B(a,a) strictly contains all Oj. The radius p(A) of the SIC of A is 
therefore strictly smaller than a. Hence, by Theorem 2.8, sin d(A, S njm ) = 
cos p{A) > cos a = (m + 1) _1 £, as claimed. □ 

The next proposition on the transition from the feasible to the infeasible 
case is similar as [15, Lemma 2.3.14]. 

Proposition 3.6 Let A = (a\, . . . ,a n ) G T n ^ m and K := — sconv(A). If 
b G K, then (A, b) := (ai, . . . , a n , b) is infeasible or ill-posed and we have 

tf{A,b) smd(b,dK) < 10tf(A). 
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Proof. W.l.o.g. A is strictly feasible. The set of solutions 

C := {x G S m | (oi,x) < 0, . . .,{a n ,x) < 0} 

is the dual of sconv(yl). Hence C = sconv( J 4). This means that a G K iff 
(a, x) > for all x G C. Therefore, we have for all a G S m , 

a K -4=>- 3x £ C (a, x) < -4=>- (ai, . . . , a n , a) is strictly feasible, (6) 

where the second equivalence follows from the assumption that C has non- 
empty interior. A similar argument shows that (ai, . . . , a n , a) is ill-posed iff 
a G dK. Therefore, we have 

d(b, dK) = min{(i(6, a) | a G S"" such that (ai, . . . , a n , a) G S n+ i )m }. 

For proving the proposition we can assume without loss of generality 
that b G K \ dK. Then (A, b) is not strictly feasible by (6). Moreover, 
since b G" dK, (A, b) is not ill-posed. Hence (A, b) is infeasible. We put now 
co := sind(b,dK) and claim that 

oj < min(6, x). (7) 

xeC 

In order to show this, suppose q G C. The equivalence (6) and b G K 
imply that cos 9 := (b, q) > 0. W.l.o.g. we may assume that — c/cos#|| 2 = 
1 — cos 2 is positive (otherwise 6 = 0, b = q, and (b, q) = 1 > to). It therefore 
makes sense to define b' := (b — qcos6)/\\b — qcos9\\. Then b' G S m and 
(V , = 0. Note that d(b, b') = ir/2 — 9. Therefore (ai, . . . , a n , 6') is feasible. 
It is either strictly feasible, in which case b' K, or ill-posed, in which case 
b' G dK (use (6)). Since b G if we conclude that d(b,dK) < d{b,b') = 
tt/2 — 8. This implies 

a; = sin d(b, dK) < cos = (b, q) 

and hence the claimed inequality (7). Moreover note that d(b,dK) < tt/2 
and to > as b G" 9if. 

Suppose now that B(p,p) is the SIC for A. Since we assume A to be 
strictly feasible t := cos p is positive. By the characterization of the GCC 
condition number in Theorem 2.8 we have t = sinci(^4, S n)jn ) = < if(A)~ 1 . 

Put ip := arcsin(jQtto) . For proving the proposition, it is enough to show 
the implication 

V(A', b') G (S m ) n+1 d{(A', b'), (A, b))<^^ (A', b') infeasible. (8) 
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Indeed, this implies (using d((A, b), £n+i,m) < tt/2, cf. Theorem 2.8) 
V(A, by 1 = sin d((A, 6), S n+ i, m ) > sin 99 = — t w = — ^(A) -1 ^, 

as claimed in the proposition. 

We argue by contradiction. Suppose there is a feasible (A', b') having 
distance at most <p from (A, b). Then there exists x' G S m such that 

K,*') < 0,...,(a^,x') < 0, (b',x'} < 0. 

Taking into account that d(a' i ,ai) < (p, we see that d(ai,x') > tt/2 — (p and 
hence (a^x') < sini^. 

We put now x := x' — Ap with A := i -1 siny?. As {a,i,p} > t, we have for 
i G [n] 

(oj, x) = (dj, x') — X(ai,p) < simp — Xt = 0. 

Note that i ^0 (otherwise t = sint/?, which is impossible). Therefore, x/||x|| 
is well-defined and lies in C. Inequality (7) implies that u>||x|| < (6, x). 

Put Ab := b'—b. Then ||A6|| < 2sm(</>/2) by our assumption d(b', b) < ip. 
We obtain 

(b,x) = (b' - Ab,x' - Xp) = (b',x') - (Ab,x') - (b',Xp) + (Ab,Xp) 
< 0+||A6|| + A+||A6||A. 

To arrive at a contradiction it is enough to verify that 

||A6|| + A + ||A6||A < uj\\x\\. 

Note that ||x|| > 1 — A, ||A6|| < 2, and to < 1. It is therefore sufficient to 
check that 

||A6|| + A + 2A < u- A, 

that is, 

||A6|| +4A < u. 

Using sin <p = 2 sin(</?/2) cos(<y?/2) we get A = t~ l simp < 2t~ 1 sin(</?/2). It is 
therefore sufficient to show that 

2 sin - + St" 1 sin - < u, 
2 2 

which is equivalent to 

(t + 4) sin I < -tu>. 

As t < 1, it is enough to show that 5 sin ^ < \toj. This is true, since by our 
assumption sin ^ < sin ip = tuj. □ 
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3.3 Feasible instances 

We provide here the proof of the part of Theorem 3.1 dealing with feasible 
instances. That is, we wish to show the claimed bound (F). 

Let A G (S m ) n , < a < vr/2, a = sina, and t > 13m(m + 1)(2<7<5 C ) _1 . 
Put e := (m+l)t _1 and ip := arcsine. We suppose that A G (S m ) n is chosen 
at random according to a distribution \i A on (S m ) n as defined in Section 2.4. 
Using Proposition 3.5 and the notation introduced there, we have 

n 

Prob{A G JT° m and V(A) > t} < Prob{A G JF° m an d ^ G T (&FQ, 

i=i 

We first bound the probability on the right-hand side for i = n by expressing 
it as an integral over A' := (a\, . . . , a n _i) of probabilities conditioned on A'. 
Note that fx A = fJ-A' x where A' := (oi, . . . , a n -i)- Moroever, A G ^ m 
iff A' G i 7 ^! m and a n G" K n , where we set now K A > := K n = — sconv(^4'), 
see (6). This implies 

Prob{A G F° n , m and a n G T (dK n ,<p)} 

Ma 



Prob{A' G m and a n G T {dK A ,,ip)} 



Ma 



= / Prob{a n G T (dK A > ,ip)} dfi AI . 

J Pan 



(9) 



A fc ^7i-l,n 



We fix A' G m and consider the properly convex set = K A r in S* m . 

The bound in Corollary 3.4 on the outer neighborhood of dK n yields 

d ur rrr(ais \\ vo\(T (dK n ,cp)nB(a n ,a)) 13m . 

Prob{a n G r o (^ n ,^)} = — — < — — sm^, 

v vo\n(a n ,a) 2a 

where v denotes the uniform distribution on B(a n ,a). The reader should 
note that e < 2cr<5 c /(13m) < a /{2m) by assumption. (We win a factor of 
two by considering B(a n ,a) instead of B(±a n ,a).) Hence, using sin p = 
e = (m + l)t -1 , we conclude 

^ , r m \ i 13m 13m(m + l) 

Prob{a n G T (dK n , p)} < — — smp - 



2a T 2at 

We assume that 13m(m + \){2at)~ l < <5 C . Hence we can apply Lemma 2.9, 
which yields 

13m(m + 1) 



Prob{a n G T (dK n , p)} < 
Pan \ 2at 
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This bound holds for any A' G J-'n-i m- We therefore obtain from (9) 
Prob{^ G T° m and a n G T (dK n ,<p)} 

The same upper bound holds for any Altogether, we obtain 

Prob{A G J-; m and Y{A) > t} < n f^^lllj ^ 
which proves Claim (F), since Probj^ G S ra ,m} = 0. 
3.4 Infeasible instances 

We start with a general observation. For A = (oi,...,a n ) G (S 1 " 1 )™ and 
1 < k < n we will write A/. := (ai, . . . , a^) and := (ai, . . . , a^). 

Lemma 3.7 Let A G (S 1 " 1 )"-, k < n, such that A k+ i be infeasible. Then 

Proof. As A k+ \ is infeasible, A must be infeasible as well. Let A' = 
(a'i,...,a' n ) be feasible such that d(A,A') = d(A,Yj n!m ) < 7r/2. Then 
A' k = (a' l5 . . . , a' k+1 ) is feasible and d(A k ,A' k ) < d(A, A'). Hence we have 
d(A k+1 ,T, k+hm ) < d(A,Z n , m ) and 

^(Ak+r)' 1 = sind(i4 fc+ i,S fc+ i )m ) < smd(A,X n , m ) = tf(A)-\ 

which was to be shown. □ 
We also need the following probabilistic lemma. 

Lemma 3.8 Let U and V be random variables taking positive values and 
xu > ol > 0, xy > P > 0, and c > 0. We assume 

Prob{[7 > x} < a ■ x~ c for x > x\j 
Prob{y>x|f/} < l3-x~ c forx>x v . 

Then we have 

Prob{[/V > x} < ca(3x~ c lnmax{ ,1} + min{a x v , f3 xfj} x~ c . 
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Proof. [23, Lemma C.l] with the functions f,g defined as 



. 1 if x < x v f 1 if x < £y 

' "^-x- c ifx>xc/ y{ ' \ l3-x~ c iix>x v 



yields 



roo 

Prob{UV>x}< f(-)(-g'( s ))ds. (11) 
Jo s 

If x > xjjxy we estimate this by 

rx/xu roo 

Prob{UV>x} < / ax~ c s c c/3s- c - 1 ds+ cps'^ds 

J Xv Jx/xrr 



' Xv Jx/xjj 
X 

caf3x~ c In ( ) + ^xfrx' 

X\jXy' 



-c 



If x < xjjxy one argues similarly. 

Finally note that (11) implies Prob{UV > x} < f °° g(j )(-/'(s)) ds, 
using integration by parts. Estimating this as before, with the roles of / 
and g exchanged, completes the proof. □ 

We provide now the proof of the part of Theorem 3.1 dealing with in- 
feasible instances, i.e., of the claimed bound (I). Fix A G S m , < a < ir/2, 
a = sin a, and t > 1. Assume A = ( 04, . . . , ct n ^j to be chosen at random ac- 
cording to (i^. Then A m+ \ is always feasible. Hence, if A = A n is infeasible, 
then there exists a smallest index k > m such that A k is feasible and A k+ \ 
is infeasible. If we denote by £ k the event 

A k feasible and A k+ i infeasible and ^(A^i) > t, (12) 

and take into account Lemma 3.7, we obtain 

n-l 

Prob{A G X„ im and V{A) > t} < ^ Prob£ fc . (13) 

fc=m+l 

For bounding the probability of £ k , a change of notation is convenient. 
We fix k and write from now on 

A := (ai,...,a k ), K A ■= -sconv{ai, . . . , a k }, b := a k+1 , 

and similarly A := (ai, . . . , d k ), b := d k+ ±. We note that A and b are chosen 
independently and at random according to (i^ and m, respectively. 
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Proposition 3.6 implies that 

( ffiA) t 

Prob£ k < Prob <^ A G T k , m and b e K A and . ,. V , ; > — 
y smd(b,KA) 10 

The first part of Theorem 3.1 tells us that 

Prob {A G ? k and If (A) > x) < k f 13m (™ + 1 ) \ x -c ^ 

provided x > xj/ := 13m(m + 1)/(2<7<5 C ). For a fixed strictly feasible A, the 
set is properly convex in S m . The bound in Corollary 3.4 on the inner 
neighborhood of dK A yields for any A G m 

r 1 I 1 1 3rr? 1 

Prob{& G K A and . — ^ > x \ A\ < , (15) 

v I sind(6, OK A ) I J 2cr x 

provided x > 2m/ a (again v denotes the uniform distribution on B(b, a)). 
Applying Lemma 2.9 yields 

p lf{" 6 K < a " d nn^SicZ) H 4 £ (M x ^ <16) 

provided x > xy := 13m/(2<7<5 c ). 

Let 1m denote the indicator function of a set M. We combine the above 
two probability estimates in (14) and (16) with Lemma 3.8, setting 

VIA) :- lr,JAmA), V(AM :_ 1 kM . 

Note that 

Prob£ fe < Prob{£/(,4) • V(A, 6) > t/10}. 
We have for t > 1 and x = t/10 



max 



fx) f 4<r 2 <^ I 

< , 1 > = max < tt, -, 1 > < t. 

\xux v J \l690m 2 (m + l)' )~ 



Let a := k (l2E!^±l)) C and £ := (i^) c . Note that a/3 = • ( 169m J m+1) 
and minjaxy, /3x^} = ^ 169 ^(™+ 1 ) ^ _ Lemma 3.8 implies that 

„ ,„ , /1690m 2 (m + l)\ c /1690m 2 (m + 1)\ c r 

Prob£ fc <c*( ^ r*lnt+( > ) r«. (17) 

Plugging in this bound into (13) finishes the proof of Theorem 3.1. □ 
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3.5 Proof of Theorem 1.2 

We will now give the proof of Theorem 1.2 by deriving estimates of the ex- 
pectation of ln'rf from the tail bounds given in Theorem 3.1. Let A G (S m ) n 
be chosen at random according to an adversarial distribution jjl^. Combin- 
ing (F) and (I) from Theorem 3.1 we have for t > 13m(m + l)(2cr<5 c ) _1 , 

Prob{tf(A) > t} = Prob{A e T n , m and V(A) > t} 

+ Prob{A G X„ im and tf(A) > t} 

< n ( ^m^+n y t-c+n ( 169 ^ ( ™ +1) ) c t- 

W ( I690^m+l) ^ t - c 

Defining D 1 := n (^^±11)' , £ 2 := £>i • (f^, Z> 3 := on A • (±f^) c 
we can write shortly 

Prob{tf(A)>t} < D 1 t~ c + D 2 t' c + D 3 t- C lnt. 

Note that D 2 > D\ and D% > D\. It is convenient to use tail estimates of 
ln^ instead of so we reformulate 

Prob{ln <g(A) > s} < D x e~ cs + D 2 e~ cs + D 3 e~ cs s, (18) 

which holds for s > In ^ 13 ™(^ +1 ) ^ _ define Sj := lnZ^^ + lntf" 1 , i.e. <5 C = 

D] /c • e- s % for i = 1,2,3. Note that s 2 = si + In > si, s 3 = 

si + c" 1 ln(cn) + In > s 1} and 

Sl = I Inn + In (±^±±>) + > In (^^) . 

The definition of <5 C in Lemma 2.9 and the inequality 
( 2 — \ -1 / 2 / — 5^ — 

( 1_ (™) m j - y ln(™/2) ' ( 19 ) 

which is shown by a small computation, lead to the following estimate of s\ 
51 ^ ln < C + ln T + ^l-(4S^)) 
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We distinguish the cases S3 > S2 and S3 < s 2 - For S3 > S2 we get with (18) 

poo 

E(ln^) = / Prob{ln V(A) > s} ds 
Jo 

rS3 roo 

< / lds+ (D 1 e- cs + D 2 e~ cs + D 3 e- cs s) ds 

JO Js a 

/•oo poo poo 

< s 3 + / D ie - cs ds+ / D 2 e- cs ds+ / D 3 e- cs sds 

J Si J S2 J SZ 

= 8$ + C' 1 81 + C' 1 81 + C~ 2 (1+083)81 

< (I + c- 1 ) s 3 + 2c- 1 + c~ 2 

= (1 + c- 1 ) si + (1 + c- 1 )^ 1 ln(cn) + (1 + cT 1 ) In 
+2C" 1 + c- 2 . 

With the estimate of si in (20) a small computation yields 

E(lnif) < (l + C - 1 )(2c~ 1 lnn + 31nm + ln(m+l)+21ni 

+(2c)" 1 In ( - 2R2 ™ ) + c- 1 + c' 1 In c + 6.5) + c" 1 . 
\m(7rm/2) / / 

This yields E(ln^) = C(c _2 ln (^)). In the case S3 < s 2 a similar argu- 
ment holds. This shows the first statement of Theorem 1.2. Finally, tracing 
the constants in the case j3 = 0, yields the asserted explicit bound. □ 

3.6 Average-case analysis 

We show here that for the uniform distribution on S m , the probability tail 
estimates in Theorem 3.1 on ^(A) can be significantly improved by essen- 
tially the same method. 

Proposition 3.9 For A G (S m ) n chosen uniformly at random we have 

Prob{^(^l)>t} = 0((m + l) 5 j Int). 
Moreover, E(\n^'(A)) = 0(\nm) as stated in Corollary 1.3. 
Proof (Sketch). A result due to Wendel [31] states that for k > m 

^' m):= (v T^F = ^ 3T ^V i )• (21) 



i=0 
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We refer to the proof in §3.3 for the uniform distribution on S m instead 
of Ha (think of a = 1). We write k instead of n. We do not need to use 
Lemma 2.9. Moreover we do not bound p(k, m) = Probj^ G J~k,m} by 1 as 
in Equation (10). Taking account of this, the proof in §3.3 shows that 

r-, , r . , ,n , 13m(m + 1) . 1 

Prob{A G F Km and Sf(A) > x) < k y — J -p(k,m) -, 

provided x > xjj := 2m(m + 1). 

We proceed now as in §3.4, using the same notation. For a fixed strictly 
feasible A = (a\, . . . , a&) we have by (15) 



Prob 



\b G K A and . >iU < 



sind(6, <9-fC0 J 2 x 



provided x > xy := 2m. Recall the definition of the event from (12). 
Similarly as for (17) we conclude with the help of Lemma 3.8 that 

Probffe < Cm 3 kp(k,m) j Int for t > e, 

where C stands for a universal constant. Using Lemma 3.10 stated below 
we get 

n— 1 4m oo 

kp(k,m) < kp{k,m)+ kp(k,m) = 0(m 2 ). 

k=m+l k=m+l fc=4m+l 

Hence, by Equation (13), 

n— 1 ^ 

Prob{,4 G Z„, m and if (A) > t} < Prob£ fc < Cm 5 - lnt 

fc=m+l 

for some constant C. It is obvious that Probj^ G J- n ,m and ^(A) > i} can 
also be bounded this way. Finally, the claimed bound on the expectation of 
ln^A) follows immediately with the help of [6, Prop. 2.4]. □ 

Lemma 3.10 We have YlV=Am+\ kp(k, m) = o(l) for m — > oo. 

Proof. Let > 4m. Wendel's result (21) implies 

, ^ . (m + 1) 2(m + l) fc m+1 

kp{k,m) < k < 



2 k 1 \ m 1 ml 2 k 
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We have k m+1 2 k < 2 k l 2 for k > Cm log to, and sufficiently large to, where 
C > is a suitable universal constant. Therefore, we get 

E, . 2(TO + f ) 1 ... , 

kp(k,m) < — — ^— ^_ = o(l)( m ->oo). 

fc>Cmlogm k=0 

The function x i— >■ x m+l 2~ x is monotonically decreasing for x > (to+1)/ In 2. 
Hence, as k > 4to, and using to! > (m/e) m we get 

f 1 (4TO) m+1 ,ev m 

— r — i— < — r - — j < 4to - . 

to! 2 fe ~ to! 2 4m ~ y 4' 

Since e/4 < 1, we conclude 

Cm log m 

kp(k,m) < 8to(to + 1) (-) m CTO log to = o(l) (to — > oo), 

fc=4m+l 

which completes the proof. □ 



4 Some spherical convex geometry 

The goal of this section is to provide the proof of Theorem 3.3, following the 
lines of [7, Theorem 1.2]. We proceed in several steps. 

4.1 Integrals of curvature and Weyl's tube formula 

For the following material from differential geometry we refer to [29]. A 
good reference for the differential geometry of convex sets is [3]. 

Let V be a smooth hypersurface in S m with unit normal vector field 
v ~. V — y S m . The principal curvatures of V at x G V are defined as the 
eigenvalues ki(x), . . . , K m -i(x) of the Weingarten map —Dv(x): T X V — > 
T X V . The ith curvature J^v,i(x) of V at x is the ith symmetric polynomial 
in the principal curvatures: 

Jty,i{x) := $ — * < 

\i\=ijei 

Interesting special cases are <%v,o{x) = 1 and 

Jfv,m-i(x) = ki(x)--- K m -i(x) = det(-Du(x)), (22) 
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which is called the Gaussian curvature of 7 at s. Let U be an open subset 
of V. In [7] the integral Hi(U) of ith curvature and the integral \/j,i\(U) of 
ith absolute curvature were defined as 

IH{U):= [ Jf V:i dV, \m\(U):= f \Jf v ,i\dV. 
Ju Ju 

Two special cases deserve special mention: no(U) = vol U equals the (m— 1)- 
dimensional volume of U. Moreover, /x m -i(V) is the integral of the Gaussian 
curvature of V. 

By a smooth convex body K in S m we will understand a convex body 
such that its boundary dK is a smooth hypersurface in S m (of type C°°) 
and its Gaussian curvature does not vanish in any point of dK. 

Let K be a smooth convex body in S m with boundary V := dK. We 
denote by v : V — > S m the unit normal vector field of the hypersurface V that 
points inwards of K. Here all the principal curvatures Kj (x) are nonnegative, 
cf. [3]. Hence the ith curvatures are nonnegative as well and therefore we 
have Hi{U) = \ni\(JJ) for any open subset U of V. 

For < ip < 7r/2 we define the (p-tube T ± (U, ip) around U as 

T-^C/, <p) := {x £ S m | 3y G U such that d(x, y) < if and 

[x,y] intersects U orthogonally at y}. 

The outer (p-tube T^-(U, <p) and inner p-tube T^-(U, p) of U are defined as 

T^U, p) := T L (U, p) \ K and Tj-{U, p) := T\U, p) n K 

In an important paper, Weyl [32] derived a formula for the volume of 
tubes around compact submanifolds of euclidean spaces or spheres. His 
result can be seen as extension of Steiner's formula on the volume of "parallel 
convex sets" in euclidean space, see also Allendoerfer [1]. Weyl's formula 
only holds for a sufficiently small radius. In [7, Prop. 3.1], it was observed 
that when replacing integrals of curvature by absolute integrals of curvature, 
one obtains an upper bound on the volume of tubes holding for any radius. 
As the above two notions of curvature coincide for boundaries of convex 
sets, we get the following result. (An inspection of the proof of [7, Prop. 3.1] 
reveals that separate bounds on the inner and outer tube hold.) 

Proposition 4.1 Let K be a smooth convex body in S m and U be an open 
subset of dK. Then we have for all < p < tt/2 

m—l 

^{vo\T^{U,p),vo\Tt(U,p)} < J2 J m,i+i(f)^(U), 

i=0 
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where J m ,i+i denotes the function defined in (4). Moreover, this upper 
bound is sharp for sufficiently small (p, cf. [32]. 

4.2 Some integral geometry 

We will need a special case of the principal kinematic formula of integral 
geometry for spheres. We denote by G the orthogonal group 0(m + 1), that 
operates on S m in the natural way, and denote by dG its volume element 
normalized such that the volume of G equals one. The following result is 
Theorem 2.7. in [7]. (For related information see [19] and [16].) 

Theorem 4.2 Let U be an open subset of a compact oriented smooth hy- 
persurface M of S m and < i < m — 1. Then we have 

Hi(U) = C(m, i) f ^(gU n S* +1 ) dG(g), 

JgeG 

where C(m,i) = (m-i-l) f" 8 : 1 ) n ° m -}? m . 

The special case i = yields an effective tool for estimating volumes, 
usually referred to as Poincare's formula: 

vol m _ 1 l/=%^/ ^(UngS^dGig), (23) 
1 Jg&G 

where C\gS l ) denotes the number of elements in U HgS 1 (note that this 
is a finite set for almost all g £ G). Here is an application of (23). Clearly, 
the given bound is sharp (consider spherical caps with radius almost tt/2). 

Corollary 4.3 Any smooth convex body K in S m satisfies vol dK < O m -\. 

Proof. Almost surely, the intersection dK n gS 1 is finite. Then it consists 
of at most two points by convexity. □ 

4.3 Integrals of curvature for boundaries of convex sets 

In this section we assume that K is a smooth convex body in S m and v is 
the unit normal vector field on dK pointing inwards of K. This means that 
for x £ V, the unit vector v{x) is uniquely characterized by the conditions 
(v, x) = and (v, y) > for all y G K. 

Lemma 4.4 We have —v(dK) = dK. 
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Proof. The characterization of u(x) implies that —v(x) G dK for x G <9if. 
For the other inclusion, let v be a unit vector satisfying —v G di'T. Then 
(f, y) > for all y £ K. Moreover, there exists x £ K such that {v,x) = 0. 
This implies x £ dK. It follows that v = v(x). □ 

The following bound is crucial for all what follows. Again, considering 
spherical caps with radius almost tt/2, shows the optimality of the bound. 

Proposition 4.5 The integral of Gaussian curvature of dK is bounded as 

fJ'm—l 

(dK) < e> m _L 

Proof. Put V := dK. By Lemma 4.4, v. V — > dK is surjective. By (22) 
we have J^v,m-i(x) = det(— Dv{x)) for x £ V. Since we assume that the 
Gaussian curvature does not vanish, the map v has no singular values. 

We claim that v is injective. Otherwise, we had v{x) = v(y) for distinct 
x,y £ V. Since (v(x),x) = and (v(y),y) = we had (i/(x),z) = for all 
z £ [x,y]. Hence v would be constant along this segment and therefore x 
would be a critical point, contradicting our asumption. 

We conclude that —v. V — > v(V) is a diffeomorphism onto the smooth 
hypersurface dK. The transformation theorem yields 

f x m _ 1 (V)= [ ^ v ,m-idV= [ det{-Dv) dV = vol dK. 
JV ' JV 

Corollary 4.3 implies now the assertion. □ 



Lemma 4.6 For a £ S m , < a < tt/2, a = sin a, and < i < m we have 

m(dKnB(a,a)) < ~ ^ O m _i a™"*" 1 . 

Proof. This is similar, but somewhat simpler than the proof of [7, Prop. 3.2]. 
The case i = m — 1 is already established by Proposition 4.5. Hence we as- 
sume i < m — 1. Let g £ G = 0(m + 1) be such that V := dK intersects 
gS l+1 transversally with nonempty intersection. We apply Proposition 4.5 
to the convex body K n gS t+l in the sphere gS t+1 , which has the smooth 
boundary V P\gS i+1 . Hence [J,i(V D gS t+1 ) < Oi. The kinematic formula of 
Theorem 4.2 applied to the open subset U := V D \nt(B(a, a)) of V yields 

m{U) = C(m, i) [ fnigU n S t+1 ) dG(g) 

Jg&G 

< C(m, i) Oi Prob{gU n S i+l + 0}. 

geG 
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Using gU C B(ga, a), this probability may be estimated as follows 
Prob{gU D S i+1 ^ $} < Prob{B(ga, a) n S i+1 + 0} 
= Prob {S(a', a) n / 0} = C^volT^ 1 , a). 

Lemma 2.1 in [7] implies vol T(S l+1 , a) = Oi + \O m -i-2 Jm,m-i-i(®) ■ Mo- 
roever, Lemma 2.2 in [7] says that 

k 

J m ,k{a) < y for 1 < fc < m. (24) 

By combining these estimates and plugging in the formula for C(m,i) from 
Theorem 4.2, the resulting expression considerably simplifies and we get 
Hi{U) < ( m ~ 1 )O m _i a" 1 -*- 1 as claimed. □ 



4.4 Proof of Theorem 3.3 

We can finally provide the proof of Theorem 3.3. We assume first that K is 
a smooth convex body in S m . Let a £ S m , < a, <p < vr/2, put a = sin a, 
e = sin ip, and let U = dK n B(a,a). By combining Proposition 4.1 with 
Lemma 4.6 we get 

m 1 / — 1\ 

volT^(C^) < ^ m . JO^-ia" 1 -*- 1 ^!^). 
i=0 ^ i ' 

Using the estimate (24) we obtain after a short calculation (put k = i + 1, 
use (T-i) = m (T ) an< ^ cons ider separately the term for k = m) 

vo\Tj-{dKnB(a,a),<p) < ^2=1 £ ( ™ j e * a™"* + -O m e m . (25) 

171 k=i ^ ' 

The same upper bound holds for the volume of T^(dK n B(a, a),ip). 
We claim that 

T (dK, ip) n B(±o, a) C T^(dK n B(±a, /?), (26) 

where /3 = arcsinminjl, a + e}. Indeed, suppose x € T o (0K, 99) n 5(±a, a) 
and let y G <9-?T be a closest point to x. Then d(x, y) < ip and [x, y] in- 
tersects orthogonally (as dK is smooth without boundary). The tri- 
angle inequality for projective distance (cf. §2.2) implies that sind(a, y) < 
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smd(a,x) + sin d(x,y) < a + e. Hence sind(a, y) < sin/3 and therefore 
y G B(±a,f3) which shows the claim. 
By combining (26) with (25) we get 

vo\(T (dK,<p)nB(±a,a)) < _^zl ^ (\ £ k {a + e) m-k + 0m£ m 

m k=i ^ ' 

We have vo\B(±a,a) > 2C m _i^, cf. [7, Lemmas 2.1-2.2]. Using this, we 
obtain 

vo\(T (dK, <p)r\B(±a,a)) fm\ / eyn-k , £ \k mO m / e\m 

vo\B(±a,a) ~ ^ \k J \ a) W 2C m _i \a) ' 

k=l 

(27) 

This shows the assertion of Theorem 3.3 for the outer neighborhood in the 
case where K is a smooth convex body. The bound for the inner neighbor- 
hood is shown similarly. 

The general case where K is any properly convex set in S m will follow 
by a pertubation argument. We define the Hausdorff distance d(K, K') of 
two convex sets K and K' in S m as the infimum of the real numbers 5 > 
satisfying K C T(K', 5) and K' C T(K, 5). This defines a metric and allows 
to speak about the convergence of convex sets. (For compact convex sets in 
euclidean space the Hausdorff distance is a well known notion.) 

Lemma 4.7 Any properly convex set K in S m is the limit of a sequence of 
smooth convex bodies. 

Proof. The euclidean version of the claim is a well known result due to 
Minkowski, see [3, §6] (or [24]) for more information. 

A properly convex set K C S m is contained in an open halfspace. For 
fixed p G S m consider now the open halfsphere S 1 ™ := {x G S m \ (x,p) > 
0} with center p and the affine space E := {x G IR m+1 | (x,p) = 1}. 
The "perspective map" it: S™ 4 4 (p,x)~ 1 x maps an intersection 

of a linear space with S m to an affine linear subspace of E and vice versa. 
Moreover, 7r maps convex sets to convex sets and vice versa. Moreover, one 
sees that it induces a homeomorphism between the set of convex subsets of 
ST and the set of compact convex subsets of E. It is easily checked that 
if K C E smooth compact convex has nonvanishing Gaussian curvature on 
the boundary, then this also holds for ir(K). The assertion follows from the 
euclidean version of our claim. □ 
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To finish the proof of Theorem 3.3 let now K C S m be a properly 
convex set and 5 > 0. By Lemma 4.7 there exists a smooth convex body K' 
such that K and K' have Hausdorff distance at most 5, which means that 
K C T(K', (5) and K' C T(if, (5). This implies K'\K C. T(dK, 5) and 

T (dif, 99) C T (dK', <p + S)U(K'\ K). 

By applying (27) to T (dK', <p+5), letting <5 0, and noting that vo\T{dK, 5) 
goes to zero, the desired assertion follows. For the inner neighborhood one 
argues similarly. □ 
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